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ABSTRACT 

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre 
(submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-lR/submm wavelengths 
over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC 
catalogue to /iM5-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IlFSCz), so that we can measure the spectral energy 
distributions (SEDs) of these objects from 60 to 850 yum. This produces a list of 17 17 galaxies with reliable associations between Planck and IRAS, 
from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no 
evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. 
We find evidence for colder dust than has previously been found in external galaxies, with T < 20 K. Such cold temperatures are found using 
both the standard single temperature dust model with variable emissivity /?, or a two dust temperature model with fS fixed at 2. We also compare 
our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find 
that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using 
more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parametric fitting. We 
thus conclude that cold (T < 20 K) dust is a significant and largely unexplored component of many nearby galaxies. 

^ Key words. Infrared: galaxies - Submillimetre: galaxies - Galaxies: ISM 



1 . Introduction tion (SED) of dust emission, combined with limited wavelength 

coverage, means that IRAS was insensitive to dust below a tem- 
Dust IS an important constituent of the interstellar medium (ISM) pgrature of ~ 30 K. Observations of dust in our own Galaxy by 
of galaxies. Whilst some properties of dust in our own and ^^e FIRAS instrument on COBE (Reach et al. 1995) found ev- 
very nearby galaxies can be studied through its absorption of ^^^^^^ ^^^^ ^^^^^^j different temperatures. This included 
starlight. It was the IRAS satellite that first allowed dust emis- ^ widespread component at 16-21 K, and another at 10-14K 
sion to be directly observed in large samples of external galax- associated with molecular clouds in the inner Galaxy. A third 
les (e.g., Devereux & Young 1990). The all-sky IRAS survey widespread colder component, at 4-7 K, was later identified with 
at 12, 25, 60 and 100 ^m in wavelength provided much new ^^e Cosmic Infrared Background (CIB; Puget et al. 1996; Fixsen 
information on the properties of dust and how this relates to ^gg^y ^one of these components would be detectable in 
other aspects of galaxies and galaxy evolution. However, the g^^^^^^j g^j^^jg^ jj^^ COBE-DIRBE observations also de- 
strong temperature dependence of the spectral energy distnbu- ^^^^^^ g^jg^^^j galaxies, finding an average dust temperature 
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of 27.6 K (Odenwald et al. 1996). 
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Observations at longer far-infrared (FIR) or submillimetre 
(submm) wavelengths from the ground (e.g., Dunne et al. 2000), 
from space (e.g.. Dale et al. 2005) or in combination (e.g., 
Willmer et al. 2009), have provided hints that cooler dust plays 
a significant role in nearby galaxies. Observations with Herschel 
of pre-selected objects (e.g., Boselli et al. 2010) or of relatively 
small fields (e.g., H-ATLAS Bales et al. 2010, covering up to 
550 deg^) provide valuable data at 250 to 500 //m, which con- 
strain the long wavelength dust properties for specific popula- 
tions. However, the availability of the Planck^ Early Release 
Compact Source Catalogue (ERCSC) provides a long wave- 
length counterpart to IRAS, allowing us an unbiased view of 
the FIR-to-submm SEDs of a large, unbiased sample of nearby 
{z < 0.25) galaxies. We are now, for the first time, able to exam- 
ine the role of cold dust for a wide range of objects in the local 
Universe. 

The discovery of the CIB (Puget et al. 1996; Fixsen et al. 
1998) has added to the importance of our understanding of 
dust in galaxies. The CIB demonstrates that roughly 50% of 
all energy generated in the history of the Universe was ab- 
sorbed by dust and reprocessed into the FIR/submm (Gispert 
et al. 2000). Deep surveys at 850yum with SCUBA (e.g.. Small 
et al. 1997; Hughes et al. 1998; Eales et al. 2000; Coppin et al. 
2006) and at nearby wavelengths with other instruments (e.g., 
MAMBO, AzTEC and LABOCA) have revealed much higher 
number counts than would be predicted by a non-evolving ex- 
trapolation of the local population. There must thus be very 
rapid evolution of the FIR/submm galaxy population, something 
confirmed by observations with ISO, (e.g.. Dole et al. 2001) 
Spitzer (e.g., Frayer et al. 2006; Bethermin et al. 2010; Clements 
et al. 2010a) and BLAST (DevUn et al. 2009). Herschel ob- 
servations have now confirmed this rapid evolution through a 
combination of number count (Clements et al. 2010c; Oliver 
et al. 2010) and luminosity function (Dye et al. 2010) studies. 
However, detailed interpretation of these results is hampered by 
our poor knowledge of galaxy SEDs in the 100-1000 fim range. 
This is demonstrated, for example, by the apparent separation in 
the temperature-luminosity plane of local //MS'-selected galaxies 
and high redshift SCUBA-selected SMGs (e.g., Clements et al. 
2010b). The origin of this separation is unclear It might rep- 
resent a genuine change in dust temperature with redshift, and 
selection biases may be partly involved, but it could also reflect 
our ignorance of the full FIR/submm SED of local galaxies. By 
properly establishing a zero redshift baseline for the dust SEDs 
of typical galaxies, the Planck ERCSC will allow the origins of 
the CIB and the nature of the galaxies that contribute to it to be 
much better determined. 

The central goals of this paper are thus twofold: to exam- 
ine the properties of a large sample of local (z < 0.25) galaxies 
to establish the range of dust temperatures and other properties 
found locally; and thus to set the local baseline against which 
higher redshift studies, and especially studies of the SMGs re- 
sponsible for the CIB, can be compared. 

The rest of this paper is organised as follows. In Section 2 
we give details of Planck's observations of local galaxies and 
the ERCSC. We also discuss the results of matching ERCSC 
galaxies to sources observed by IRAS. In Section 3 we present a 



' Planck (http://www.esa.int/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries 
France and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



comparison of the ERCSC with existing data from ground-based 
observatories. Section 4 presents the results of fitting parametric 
models to the dust SEDs of ERCSC galaxies, while Section 5 
discusses the results of physical template fitting. Finally, we 
draw conclusions in Section 6. Throughout this paper we as- 
sume a concordance cosmology, with //q = 70kms"'Mpc"', 
Qa = 0.7 and Qm = 0.3. 

2. Planck Observations of Nearby Galaxies 

2. 1 . The Planck Mission 

Planck (Tauber et al. 2010; Planck Collaboration 2011a) is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. The Low Frequency 
Instrument LFI; (Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 201 1) covers the 28.5, 44.1, and 70.3 GHz bands, 
with amplifiers cooled to 20 K. The High Frequency Instrument 
(HFI; Lamarre et al. 2010; Planck HFI Core Team 2011a) cov- 
ers the 100, 143, 217, 353, 545, and 857 GHz bands, with 
bolometers cooled to 0.1 K. Polarization is measured in all but 
the highest two bands (Leahy et al. 2010; Rosset et al. 2010). 
A combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics (Planck Collaboration 201 lb). Two data processing centers 
(DPCs) check and calibrate the data and make maps of the sky 
(Planck HFI Core Team 2011b; Zacchei et al. 2011). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration, 201 Ih-z. 

2.2. The Planck Early Release Compact Source Catalogue 

The Planck ERCSC (Planck Collaboration 201 Ic) provides po- 
sitions and flux densities of compact sources found in each of the 
nine Planck frequency maps. The flux densities are calculated 
using aperture photometry, with careful modelling of Planck's 
elliptical beams. The colour corrections for sources with spec- 
tral index a - -0.5 (using the convention Sy oc v") are 1.017, 
1.021 and 1.030, respectively, for the 28.5, 44.1, and 70.3 GHz 
LFI channels. Flux densities taken from the ERCSC should be 
divided by the appropriate colour correction to give the correct 
flux values for an assumed narrow band measurement at the 
central frequency. For frequencies from 28.5 to 143 GHz, com- 
pact sources have been detected using a version of the "Powell 
Snakes" techniques (CarvaUio et al. 2009); for details see Planck 
Collaboration (2011c). In the four higher frequency channels, 
sources were located using the SExtractor package (Bertin & 
Arnouts 1996). Sources detected in one or more of the frequency 
maps were then put through a further set of secondary selection 
criteria; these are discussed in detail in Planck Collaboration 
(2011c). The primary criterion utilized was a Monte Carlo as- 
sessment designed to ensure that >90% of the sources in the 
catalogue have a flux accuracy of at least 30%. 

2.3. Matching The ERCSC to IRAS Data 

To understand the FIR Spectral SEDs we need a combination 
of data at long wavelengths, provided by the Planck ERCSC, 
and data near the peak of a typical galaxy dust SED at about 
lOOjum. The best source for the latter information is the IRAS 
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Fig. 1. Histogram of offsets between ERCSC and IIFSCz posi- 
tions. 



test the validity of the association with the IRAS source. This in- 
cluded inspection of the Sky Survey postage stamps provided in 
NED. Associations were accepted as real if the associated galaxy 
had a blue {g or B) magnitude brighter than 16. The surface- 
density of such galaxies leads to the probability of a chance as- 
sociation being ~3%. For sources where there was both an IRAS 
and a 2MASS association, this limit was relaxed to B = 17 (or 
K ~ U). Of the 88 ERCSC-IIFSCz associations with positional 
offsets 3-5' (and with spectroscopic or photometric redshifts) 
72 were associated with bright galaxies, two were associated 
with a second FSC source having cirrus-like colours and are pre- 
sumed to be cirrus, and the remaining 14 are classified as pos- 
sible galaxy associations (these are excluded from the reliable 
galaxy catalogue used here for further analysis). 

The second category of ERCSC-IIFSCz associations which 
we scrutinised in detail were those for which there is no redshift 
in the IIFSCz. There were 165 of these and the NED associa- 
tions suggest that 38 are bright galaxies, seven are cirrus, two are 
bright planetary nebulae, and the remaining 118 are classified as 
possible galaxy associations (and excluded here). We were left 
with 1717 reliable ERCSC-IIFSCz galaxy associations of which 
337 are flagged as extended in the ERCSC. 1597 of these 1717 
objects have spectroscopic redshifts. 



all-sky FIR survey, and the most recent analysis of the IRAS 
Faint Source Catalogue (FSC) is provided by Wang & Rowan- 
Robinson (2009) in the Imperial IRAS FSC redshift survey 
(IIFSCz). This was constructed using IRAS FSC sources, all 
of which are at \b\ > 20°, with IRAS colours used to exclude 
stars and cirrus sources (with 5(100)/5(60) > 8). The NASA 
Extragalactic Database (NED) was then used to find spectro- 
scopic redshifts for the resulting FSC source list, and to asso- 
ciate the sources with SDSS (where available; York et al. 2000) 
and 2MASS (Skrutskie et al. 2006) galaxies to find photometry 
at 0.36-2.2 fim. This photometry was then used to estimate pho- 
tometric redshifts for sources without spectroscopic redshifts. 

The starting point for matching the ERCSC to the IIFSCz is 
the 9042 sources detected at 857 GHz by Planck. This is then 
restricted to the 5773 sources at \b\ > 20° for which there will 
be FSC data. Associations of ERCSC sources with IIFSCz were 
looked for using a search radius of 5'. The histogram of posi- 
tional offsets is shown in Fig. 1. The bulk of associations have 
offsets within 2'. Even at 5' there is no steep increase in the num- 
ber of associations which would be indicative of a large fraction 
of chance associations. On the basis of source surface-density, 
the chance of a random association with an IIFSCz source within 
3' is 1.6%, and within 5' it is 4.5%. A total of 1966 associations 
were found within 5'. There were 106 cases where an ERCSC 
source picked up an association with more than one IIFSCz 
source. We examined these cases carefully to ensure that only a 
single association was accepted. Generally the nearer association 
was preferred. Where the offsets of the two associations were 
comparable, the brighter IRAS source was preferred. There are 
20 cases where two bright galaxies less than 5' apart have gener- 
ated a single ERCSC source, for which there may be a significant 
contribution from both galaxies to the ERCSC flux. These con- 
fused sources would benefit from additional observations with 
ground-based submm instruments to determine the contribution 
of each component to the submm emission detected by Planck. 

The remaining ERCSC-IIFSCz associations were further 
scrutinised as follows. Firstly, for sources with positional off- 
sets between the two catalogues in the range 3-5' all NED as- 
sociations within 5' of the ERCSC positions were examined to 



2.4. ERCSC Sources not associated with IIFSCz Sources 

Fig. 2 shows the sky distribution of ERCSC sources at \b\ > 20°, 
with sources flagged as extended in the ERCSC shown as blue 
filled hexagons, and point-sources shown in black. Associations 
with the IIFSCz are shown as red circles. The extended sources 
not associated with IIFSCz sources have a strikingly clustered 
distribution, which matches the areas of our Galaxy with strong 
cirrus emission, as evidenced by IRAS lOOyum maps and by the 
ERCSC cirrus flag (values > 0.25). We presume the majority of 
these are cirrus sources and not extragalactic. 

To test this further, we looked for NED associations with 
all 444 extended ERCSC sources lacking IIFSCz associations at 
\b\ > 60°. Only 12 were found to have associations with bright 
(b, g < 16) galaxies. Extrapolating to \b\ - 20-60°, we estimate 
that a further ~ 50 of these extended non-FSC sources will be 
bright galaxies. The remainder of the 3431 extended non-FSC 
sources at \b\ > 20° are presumed to be Galactic cirrus. 

A ridge of non-FSC point sources can be seen in Fig. 2 
at b ~ 70°, / ~ 120-230°. These correspond to one of the 
IRAS coverage gaps. We examined NED associations for all 482 
ERCSC point-sources not associated with IIFSCz sources. 32 
were found to be associated with Local Group galaxies (M31, 
SMC and WLM, with 28 in the LMC), 123 are bright galaxies, 
27 are associated with IRAS FSC or PSC Galactic cirrus sources, 
and 10 are bright stars or planetary nebulae. Most of the bright 
galaxies lie in the IRAS coverage gaps. The remaining 289 are 
classified as possible galaxy associations (and excluded here). 

To summarise, we have found a net total of 1884 definite 
associations with galaxies. These constitute an ERCSC galaxy 
catalogue. A further 419 sources are not associated with bright 
galaxies, but there are grounds for thinking they could be extra- 
galactic sources. Some have IIFSCz associations, but there are 
too many possible faint optical or near-IR galaxy counterparts to 
be confident which might be associated with the ERCSC source. 
Some of these 419 sources are almost certainly fainter galaxies, 
although many could turn out to be cirrus. Improved submm or 
FIR positions are needed, either via Herschel or ground-based 
telescopes, to identify these sources reliably. 
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Fig. 2. Sky plot of ERCSC sources in galactic coordinates. Black filled hexagons are ERCSC point-sources and blue filled hexagons 
are ERCSC sources flagged as extended. Red hexagons are sources associated with IIFSCz IRAS ESC galaxies, after scrutinising 
suspect categories with NED (and excluding some, as described in the text). Green hexagons are ERCSC sources not associated 
with IIFSCz, but associated with bright galaxies in NED (only for \b\ > 60° for extended sources). 



Following this identification analysis we restrict ourselves to 
those galaxies with reliable IIFSCz associations and with de- 
tections in the 857 and 545 GHz bands at significance of So- 
ot greater, as well as detections in the 353 GHz band of 3cr or 
greater This amounts to a total sample size of 595 galaxies. 

2.5. Cirrus Contamination 

Our analysis of the non-IIFSCz-identified ERCSC sources in 
Section 2.4 led us to the conclusion that sources which are 
extended in the ERCSC are a result of cirrus structure in our 
own Galaxy, or at best are a combination of cirrus structure 
with flux from a galaxy. Of the 595 reliably detected IIFSCz- 
identified ERCSC sources, 127 are listed as extended in the 
ERCSC. We test these objects for the possibility of cirrus con- 
tamination by examining the amplitude of the local cirrus fluc- 
tuations in the 100 yum cirrus maps of Schlegel et al. (1998). We 
adopt this approach since regions of greatest cirrus fluctuation 



are those most likely to cause problems for point source detec- 
tion in the ERCSC. We measure the cirrus RMS in a 3 x 3 array 
of points, separated by 0.1° and centred on the position of the 
ERCSC source. Since cirrus emission is likely to have cooler 
FIR-to-submm colours than the integrated emission of an exter- 
nal galaxy, we then look for any correlation between cirrus RMS 
and the 60 yL(m-to-857 GHz colour We plot this relation in Fig. 3. 

As can be seem from Fig. 3, there appears to be a correla- 
tion between colour and cirrus RMS for the sources classified 
as extended in the ERCSC. We conclude that the ERCSC fluxes 
for these sources are partially contaminated with cirrus emission 
from our own Galaxy. We thus exclude these 127 sources from 
further analysis. Of the remaining 468 non-extended ERCSC 
sources, fewer than 10 lie in the region of this correlation. These 
sources are retained for the following analysis, but any conclu- 
sions that come solely from these specific sources will be treated 
with caution. 
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Fig. 3. 60 //m to 857 GHz (i.e., 350//m) colour plotted against 
the cirrus RMS at 100/Ym in IRAS. ERCSC sources classified as 
point-like are shown as open diamonds, while extended sources 
are shown as solid dots. Note that the extended sources show a 
clear correlation between colour and cirrus RMS, indicating that 
these sources are likely to be contaminated by cirrus emission. 

More generally, this analysis highlights one of the issues that 
must be faced when using the ERCSC catalogue. Anyone wish- 
ing to cross-match Planck sources, especially those detected at 
high frequencies, with sources at other wavelengths, needs to 
take great care in ensuring that the ERCSC fluxes are not con- 
taminated by cirrus emission. 

3. Comparison to Existing Submm Data 

3.1. Galaxies detected with SCUBA 

The largest studies of cool dust in external galaxies to date have 
been associated with the SCUBA Local Universe Galaxy Survey 
(SLUGS) and its extensions (Dunne et al. 2000; Dunne & Eales 
2001; Vlahakis et al. 2005; Clements et al. 2010b). These en- 
compass a total of about 250 objects that were observed with 
SCUBA. The targets were selected on the basis of IRAS flux B- 
band optical magnitude or FIR luminosity. Most of the objects 
were detected only at 850 yum (i.e., not also at 450 jum), allowing, 
with the IRAS data, only a single component {T,/3) fit - where 
the SED is described as Sy cc v^B{v, T), with B{v, T) being the 
Planck function, and the parameters T and /? representing tem- 
perature and dust emissivity index, respectively. A small fraction 
of SLUGS galaxies were also detected at 450 fim, allowing for 
the existence of a second, cooler, dust component to be assessed. 
For these objects, and more recently for an ultra-luminous IR 
galaxy (ULIRG) sample, Dunne & Eales (2001) and Clements 
et al. (2010b) found some evidence for a colder dust contribu- 
tion. 

The presence of colder dust can be inferred from colour- 
colour diagrams when two submm flux densities are available. 
We show the SLUGS sources and the ERCSC sources (after 
colour corrections to the Planck flux densities and a suitable 
scaling has been applied to convert from Planck 857 GHz flux 
density to the SCUBA 450 //m band) in Fig. 4. As can be seen, 
the Planck galaxies lie on the same broad trend as the SLUGS 
galaxies (with the exception of a small number of objects dom- 
inated by a non-thermal AGN component, such as 3C273 and 
3C279). The ERCSC sources, though, extend the trend to cooler 



FIR/submm colours than were found for the SLUGS objects, 
suggesting that the galaxies detected in the ERCSC contain 
cooler dust than was detected in the majority of SLUGS sources. 

3.2. CO Contamination 

One factor that has complicated the interpretation of ground- 
based submm observations of galaxies has been the presence of 
CO emission lines within the submm passbands that make a sig- 
nificant contribution to the continuum flux. Seaquist et al. (2004) 
estimated that the CO (3-2) line contributed an average 25% of 
the flux received in the SCUBA 850 //m continuum passband 
for galaxies observed in the SLUGS survey, with the range of 
flux contributions going from 10-45% for the subset of SLUGS 
galaxies for which CO(3-2) observations were available. The 
SCUBA 850 yum filter has a bandwidth of -30 GHz. The Planck 
353 filter is significantly broader, at ~90GHz, so the line contri- 
bution will be correspondingly smaller at ~8% on average. This 
contamination fraction can be checked using observed values of 
integrated CO 3-2 line fluxes for a variety of objects from Bayet 
et al. (2006) and matching them to continuum observations of 
similar beamsize to the CO observations. We find contamination 
fractions of 2% for Arp220, 6% for Mrk23 1 and a worst case ex- 
ample in the central region of NGC 253, where a contamination 
of 1 1 % is calculated. 

We can extend this analysis to other Planck bands using 
NGC 253 as a worst case since, unlike most others, this object 
has been observed over the full range of CO transitions acces- 
sible from the ground. We find that the higher frequency bands 
have a reduced level of contamination compared to the 353 chan- 
nel, with <1% at 857 and 6% at 545 (assuming a flat spectral 
line energy distribution to estimate the contribution of the CO 
5-4 line that is inaccessible from the ground). More normal ob- 
jects than NGC 253, not dominated by an ongoing starburst, will 
have an even smaller level of CO copntamination than this. At 
lower frequencies, though, the contamination can become more 
serious. CO 2-1 could contribute as much as 21% of the con- 
tinuum flux in the 217 band in the inner regions of NGC253. 
In the 100 band for this object the CO line could contribute as 
much as 75% of the flux of the thermal continuum. NGC253 is 
a worst case scenario, so more typical sources would of course 
suffer much less contamination. However, very few galaxies are 
detected by Planck solely in thermal emission in this band, and 
the few that are detected are bright nearby objects with substan- 
tial archival data that can allow a direct assessment of the CO 
contribution. 

Our conclusion from this analysis is that the CO contribution 
to the continuum flux is likely to be smaller than other sources of 
uncertainty for generic ERCSC-detected galaxies except for the 
small number which are detected in the 217 band. Flux excesses 
detected in this band alone might thus result from CO emission 
rather than from any putative very cold dust component. 



4. Parametric IVIodels of Dust SEDs 

4.1. Fitting Metinod 

Given the Planck and IRAS flux data described in Section 2, with 
appropriate colour corrections applied to the Planck flux densi- 
ties, we model the underlying signal in observed frequency band 
V as one or more grey-body sources with flux density 

G(v; T,I3) oc yPByiT) , (1) 
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Fig. 4. Colours for the ERCSC galaxies (black dots) compared 
to those found for SLUGS galaxies (red; Dunne & Bales 2001; 
Vlahakis et al. 2005) and ULIRGs (blue; Clements et al. 2010b). 
A flux correction factor of 0.506 has been applied to the Planck 
857 GHz flux densities to extrapolate them to the SCUBA 
450/^m band. Only sources detected at > 3cr in the 353 GHz 
(850 //m) band and at > 5cr in the 857 GHz band (the require- 
ment for inclusion in our analysis) are shown. The four points 
above the general trend in the lower left of the diagram are the 
well-known non-thermal dominated sources 3C279, 0537-441, 
OJ+287 and 3C273. The SLUGS point with the lowest 60 ^im to 
SCUBA flux ratios corresponds to the galaxy IC979; it is offset 
from the general correlation for SLUGS and ERCSC galaxies, 
and Vlahakis et al. (2005) note that its IRAS flux densities should 
be treated with caution. 



where By{T) is the Planck function for blackbody flux density. 
We fit the data dv to one-component models of the form 



dy^AG[v{\ +zy,T,f3]+ny , 



(2) 



or to two-component models with a fixed j3 - 2 grey-body expo- 
nent. 



dy = AiG[v(l + z); Ti,2]+ A2G[v(l + z); Tj, 2] + ny 



(3) 



In these equations, A or A,- is an overall amplitude for each com- 
ponent, and the factor of (1 + z) converts from rest-frame fre- 
quency to observed frequency for an object at redshift z- The 
noise contribution is given by ny, which we model as a Gaussian 
with variance crj. For the Planck channels, the determination 
of the noise contribution is described in Planck Collaboration 
(201 Ic). For IRAS, the detections are classified in the IIFSCZ of 
Wang & Rowan-Robinson (2009) into (1) good detections, for 
which we take cry = Q.ldy; (2) marginal detections, for which 
we take o-y - Q.Sdy; and (3) upper limits, for which we take 
o"v to be the reported upper limit, and dy = 0. As mentioned 
above, we only consider sources with detections in the 857 and 
545 bands at significance of at least 5cr or greater and in the 353 
band of at least 3cr. 

Thus, the parameters of our model are some subset of the 
Ai, Ti,p, depending on which model we fit. We use a Bayesian 
Markov Chain Monte Carlo (MCMC) (e.g., Lewis & Bridle 
2002; Jaynes 2003) technique to probe the parameter space; with 
our Gaussian noise, this is equivalent to an exploration of the 
surface, albeit with a nonlinear parameterization. We require a 



< < 3 and 3K < T < 100 K with a uniform prior proba- 
bility between those limits (detections of very low temperatures, 
T < 10 K, are actually dominated by non-thermal emission). We 
adopt a uniform prior on In A,, as it ranges over many orders of 
magnitude for sources of widely varying absolute luminosities 
and distances. 

The MCMC engine first creates a 15,000-sample Markov 
chain, varying one parameter at a time, using this to find 
an approximately-orthogonal linear combination of parame- 
ters, with which a subsequent 100,000-sample chain is run. 
Convergence is assessed by re-running a small number of chains 
from a different starting point and checking for agreement to 
much better than one standard deviation in all parameters. 

We calculate an approximation to the Bayesian evidence, or 
model likelihood (Jaffe 1996; Jaynes 2003) in order to compare 
the two-temperature and one-temperature fits. The evidence is 
calculated as the average of the likelihood function over the prior 
distribution; we approximate the likelihood as a multivariate 
Gaussian function of the parameters centred at the maximum- 
likelihood MCMC sample with covariance given by the empir- 
ical covariance of the samples (this approximation ignores the 
prior on the amplitude of the individual grey bodies). 

In Figures 5, 6 we show sample output from our MCMC runs 
for different objects and models, along with the measured SEDs 
and fits. For objects such as F01384-7515 in Figure 5, if we in- 
stead perform a two-temperature fit, it prefers the amplitude of 
the second temperature component to be many tens of orders of 
magnitude below the first, and gives temperature values consis- 
tent with the one-temperature fit; this indicates, along with the 
approximate evidence discussed above, that a one-temperature 
model is strongly preferred. 

4.2. Results from Parametric Fits 
4.2.1. Single Component Fits 

Figure 7 shows the T plane for parametric fits to all 468 reli- 
ably identified non-extended sources within the ERCSC-IIFSCz 
cross match whose flux densities pass our S/N ratio criteria for 
inclusion, together with similar single temperature parametric 
fits from Dunne et al. (2000) and Clements et al. (2010b). As 
can be seen, the Planck ERCSC sources overlap with the SLUGS 
galaxies but extend to cooler temperatures and flatter, i.e., lower 
/3, SEDs. The median parameters for the Planck sources are T = 
26.3 K with temperatures ranging from 15 to 50 K, and fi-l.l 
compared to corresponding values from 104 SLUGS galaxies 
Dunne et al. (2000) of T = 35 K and /3=1.3 and of T = 41 K and 
y6=1.6 for 26 ULIRGs Clements et al. (2010b). This confirms the 
result from consideration of Planck-IRAS colours in Fig. 4 that 
we are seeing cooler dust in the ERCSC-IIFSCz galaxies. 

There are ten sources common to the ERCSC-IIFSCz cross- 
matched catalogue and the SLUGS studies. The fits using Planck 
data and using SLUGS data for all but two of these sources 
are in good agreement. The two exceptions are NGC7541 and 
NGC5676. The likely cause of the disagreements in these cases, 
is the presence of a close companion IRAS source to NGC7541 
(NGC7537 3.1' separation, and thus only -0.7 beam FWHM 
away), so that the Planck flux density is likely over-estimated, 
and extended IRAS emission in NGC5676, making the IRAS ESC 
flux densities used in our analysis underestimates. 

The position of our galaxies in the Luminosity-Temperature 
plane is an important question since it relates to claims of evo- 
lution in the dust properties of galaxies. It has previously been 
suggested that high redshift, high luminosity SMGs have lower 
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Fig. 5. Left panel: Samples from the SED likelihood function (y^) for FO 13 84-75 15 and the one-temperature model. The bottom row 
shows the one-dimensional (marginalized) posterior for the parameters (log A, jS, T), and the other panels show all two-dimensional 
marginal distributions. Right panel: the data points for this object, along with the one-temperature model for the maximum likelihood 
sample (yellow curve; T - 30, K, /3 - 1.7) and the model determined by the mean of the samples of each parameter (red curve; 
T = (32 + 5) K, y6 = 1 .6 + 0.4 where the errors are the variances). Note the logarithmic axes which make interpretation of the error 
bars difficult. 




log.Ai Ti log.42 wavelength [/im] 



Fig. 6. As in Figure 5, but for F00322-0840 and the two-temperature model with fixed -2. The maximum-likelihood temperatures 
are 20 K and 49 K the means and variances are (20 + 0.8) K and (44 + 8) K. 



dust temperatures and higher dust masses than more nearby ob- 
jects Yang et al. (2007). Comparison of SMGs from Chapman 
et al. (2005), Coppin et al. (2008) and Kovacs et al. (2006) with 
more local galaxies from Dunne et al. (2000) and local ULIRGs 
Clements et al. (2010b) confirms this effect. Claims have been 
made that sources selected at longer wavelengths than the 60 fim 
typical of IRAS derived samples (e.g., Symeonidis et al. 2009), 
Patel et al. (in prep)) show less of a separation between the lo- 
cal sources and the higher redshift SMGs. Much of this work is 
hampered by the poor sampling of the dust SEDs of local objects 
at wavelengths between 100 and 850 yum. Recent results from 
Herschel Amblard et al. (2010) and BLAST Dye et al. (2009) 
have begun to fill the gap between local IRAS galaxies and the 
SMGs, suggesting that our view of dust temperatures in local ob- 
jects are biased to warmer temperatures through our dependence 
on IRAS flux densities. In Fig. 8 we show the positions of Planck 
galaxies on the L-T plane. As can be seen, the gap between the 



local galaxies and the SMGs is starting to be filled by the Planck 
objects. The Planck fluxes for one of these sources may include 
some contamination from cirrus, but the rest lie in areas of nor- 
mal to low Galactic cirrus noise and should thus be fully reli- 
able. The large area coverage of Planck is particularly important 
in this as it allows us to probe generic L/ULIRG-class objects 
(LpiR > IO^Lq) rather than having to rely on pre-selected IRAS 
bright sources as in Clements et al. (2010b). Herschel observa- 
tions, which are also beginning to show the gap being filled, do 
not cover enough area to include many such L/ULRG objects in 
the local Universe Amblard et al. (2010). We find several cool 
{T < 30 K) ULIRGs that have very similar characteristics to 
SMGs. The issue of the apparent distinction between local galax- 
ies and high-z SMG population thus seems to be approaching 
resolution. 

While the temperature distribution of our objects is consis- 
tent with what has been seen elsewhere, we find that some of our 
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Fig. 7. Temperature-/? correlation for ERCSC-IIFSCz matched 
galaxies (solid dots) together with data from Dunne et al. (2000) 
and Clements et al. (2010b) for SCUBA observed sources (tri- 
angles and squares respectively). The four sources with fi-Q and 
T < lOK at the left of the diagram are the non-thermal domi- 
nated sources 3C279, 0537-441, OJ+287 and 3C273. 



galaxies have (3 more than 3cr less than 1 . One possible cause for 
this might be the 353 flux density being affected by Eddington 
bias (see e.g., Teerikorpi (2004)), leading to 'flux-boosting' of 
lower significance detections, since we accept these fluxes down 
to 3cr. This is tested by repeating the fits using only the IRAS, 
857 and 545 flux densities. While there are small differences in 
fits to individual objects resulting form the exclusion of the 353 
flux densities, the general distribution remains the same, com- 
plete with the low /? sources. We thus conclude that this is a real 
effect and not due to Eddington bias or any other issue related to 
the 353 flux densities, such as contamination by emission from 
the CO 3-2 molecular line. While such low fi values are not ex- 
pected in simple models of dust, it is suggestive that the SEDs 
can be better fit by a parameterization that uses a mixture of dust 
at two temperatures, as suggested by Dunne & Eales (2001) and 
which is a good fit to our own galaxy as seen by COBE Reach 
et al. (1995). We investigate this by applying two component fits 
to the dust SEDs. 



4.2.2. Two Temperature Fits 

We carry out two temperature component fits on our sources, as- 
suming j3=2 for both components, and present the temperature- 
temperature plot in Fig. 9. We also use the Baysian evidence 
calculated during the fitting process (Jaffe 1996; Jaynes 2003) 
to determine how many of our sources show evidence for a two 
component fit above that of the single component {T,fi) fit. We 
find that the two component fit is favoured in most cases, with 
425 objects giving a higher evidence for this model and only 43 
preferring the single component fit. Once again we test the pos- 
sibility that issues with the lower significance 353 GHz flux den- 
sities might bias these fits by repeating the analysis with these 
fluxes excluded. While this increases the uncertainties in the fits, 
as with the {T,fi) fits we find that the exclusion of the 353 GHz 
fluxes makes no systematic difference to the temperatures found. 

We find 17 galaxies fit by models containing a dust com- 
ponent with temperatures as low as 10 K. Four of these are the 



Fig. 8. The Temperature-Luminosity plane showing a variety of 
FIR populations. Open squares are SMGs from Chapman et al. 
(2005), open triangles are SMGs observed with SHARCII by 
Coppin et al. (2008) and Kovacs et al. (2006) . + signs are 
the SLUGS sources from Dunne et al. (2000), x are interme- 
diate redshift ULIRGs from Yang et al. (2007), ULIRGs from 
Clements et al. (2010b) are open circles. P/fl«c/t-ERCSC-IIFSCz 
galaxies are shown as solid dots. As can be seen the previous 
apparent distinction between the local FIR populations and the 
SMGs is weakened by sources from this work lying in the same 
region as the SMGs and by filling in some of the gap between the 
populations. The four aberrant sources with T <15K and high L 
in the bottom right of the plot are the non-thermal dominated 
sources 3C279, 0537-441, OJ+287 and 3C273. 



bright blazars and are thus dominated by non-thermal emission, 
while one is a source that might still contain some cirrus contam- 
ination. We thus find at least 13 galaxies which appear to contain 
very cold dust. Such dust has previously been found in our own 
galaxy (e.g. Reach et al. 1995) but has not been seen before in 
large scale extragalactic surveys. The small number of sources 
which have Tcoid lower than 7 K are all dominated by nonther- 
mal emission. The details of the temperature-temperature plot in 
Fig. 9, with small scatter and the temperature of the hot com- 
ponent being largely independent of that of the cold component 
up to Tcoid ~ 18 K will likely have implications concerning the 
relationship between hot and cold dust components. 



4.3. The Broader ERCSC-IIFSCz Sample 

The parametric fitting reported here only concerns the 468 non- 
extended ERCSC-IIFSCz matched galaxies which are detected 
at 5cr or greater in the 857 and 545 bands and at 3cr or greater 
at 353. This ensures that uncertainties in the fluxes do not pre- 
clude a good fit to the SED, but brings the risk that we might be 
missing a significantly different subclass of object in the remain- 
ing 1 122. We have thus applied our fitting methods to this whole 
sample, regardless of ERCSC S/N beyond the basic detection re- 
quirement of a 5cr or greater detection at 857 GHz. While there 
are larger errors bars on the fitted parameters we find no indica- 
tion that galaxies in this larger sample have a different range of 
dust properties to those discussed above. 
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use two further cirrus templates with 0=1 and 0.1, which cor- 
respond to significantly cooler dust than in the standard cirrus 
template. The starting point for the present analysis is the tem- 
plate fit for each object given in the IIFSCz Catalogue Wang 
& Rowan-Robinson (2009). This was done by fitting the op- 
tical and near-IR fluxes with an optical galaxy or QSO tem- 
plate. The IRAS data were then fitted with one of the original 
four Rowan-Robinson et al. (2008) templates. This model is then 
compared to the Planck data. In almost all cases additional com- 
ponents are needed, since the //MS'-based predictions underesti- 
mate the submm flux densities provided by Planck. More weight 
is given to 857 and 545 flux densities in this process since they 
are generally at higher signal to noise and are thus less sub- 
ject to Eddington bias effects. The 353GHz and, where detected, 
217GHz flux densities are also subject to contamination by CO 
3-2 and 2-1 emission, respectively, which may also contribute to 
an excess flux in these bands beyond what might be expected 
from the continuum fit. 



Fig. 9. The Temperature-Temperature plane for two tempera- 
ture component fits for ERCSC-IIFSCz matched galaxies, to- 
gether with data from Dunne & Bales (2001) and Clements et al. 
(2010b) for SCUBA observed sources (triangles and squares re- 
spectively). Only those ERCSC-IIFSCz sources where the two 
temperature model is preferred and with a reliably determined 
cold dust component temperature (i.e., T/o-{T) > 4) are plotted. 
The sources with the coldest T^oid < 7 K in this plot are domi- 
nated by non-thermal emission. 



5. Physical Models: Template Fitting 

The FIR and submm spectral energy distributions of galaxies 
from the IRAS, ISO and Spitzer surveys have been successfully 
modelled with a small number of templates (Rowan-Robinson 
(1992), Rowan-Robinson et al. (2005), Rowan-Robinson et al. 
(2008)). However the submm data available in such studies is 
quite limited and we expect to get a much better understanding 
of cool dust in galaxies with the data from Planck. Already with 
Herschel, there is evidence for unexpected quantities of cold dust 
in some galaxies Rowan-Robinson et al. (2010). 

The Rowan-Robinson et al. (2010). study of the SEDs of 
i/eric/ieZ-detected galaxies used SPIRE flux densities extend- 
ing from 250 to 500 jjm combined with pre-existing data from 
SWIRE at Spitzer and optical wavelengths. Through the combi- 
nation of Planck data with the IIFSCz, the galaxy sample consid- 
ered here includes data from the optical to IRAS fluxes, and then 
the Planck data has flux densities at 350, 550 and 850yum (857, 
545 and 353 GHz). Some of our objects are also detected at 1.4 
mm (217 GHz) which is included in our analysis if available. The 
range of wavelengths available with the ERCSC-IIFSCz sample 
is thus broader than that available through Herschel. This en- 
ables us to place better constraints on the role and importance of 
cold dust in these objects. Our sample is also much larger than 
the 68 objects considered in Rowan-Robinson et al. (2010), so 
we can better determine the variety and overall statistics of the 
SEDs of local galaxies. 

The templates used in fitting IRAS, ISO and Spitzer data are 
(1) a cirrus (optically thin interstellar dust) model characterised 
by a radiation intensity = I{galaxy)/I{IS RF) = 5, where 
I(ISRF) is the intensity of the radiation field in the solar neigh- 
bourhood, (2) a normal M82-like starburst, (3) a higher optical 
depth Arp220-like starburst, (4) an AGN dust torus. In fitting 
the SEDs of Herschel galaxies, Rowan-Robinson et al. (2010) 



5.1. Results from Template Fits 

In Figure 9 we analyze the SEDs of the archetypal nearby galax- 
ies, M51, MlOO, M2 and Arp 220. M51 and Ml 00 are modefled 
with two cirrus templates with (p-l (solar neighbourhood) and 5 
(Galactic Centre), and with a modest M82 starburst component. 
M82 itself needs an additional component of cool cirrus {(p-l) 
as well as the M82 template of Efstathiou and Rowan-Robinson. 
Finally Arp220 is modelled extremely well over all infrared and 
submm wavelengths by the Arp220 template used by Efstathiou 
and Rowan-Robinson. The models for these galaxies by Silva 
et al. (1998) are also shown and perform well, especially for 
M51 and Arp220. 

Fig. 10 shows fits to galaxies with detections in 16 photo- 
metric bands: 5 optical bands (SDSS), 3 near infrared bands 
(2MASS), 4 mid and FIR bands (IRAS) and 4 submm bands 
(Planck). The blue curve is the solar neighbourhood cirrus tem- 
plate {(p-l) and contributes significantly to the SEDs of 7 out 
of the 8 galaxies. Dust grains in this component are in the range 
15-20K, depending on grain radius and type (Rowan-Robinson 
1992). Dust masses are in the range 10^ - 3.10**Mo. 

To get an overview of the whole sample. Fig. 12 shows the 
distribution of 545 GHz (550 jum) flux density versus redshift for 
galaxies well detected at 350-850//m, and with spectroscopic 
redshifts. The loci of galaxies with an Arp 220 template at lu- 
minosity LiR - 10'^ and lO'^L© are shown. With the restric- 
tions to point-sources with good detections in the three highest 
frequency Planck bands, and to galaxies with spectroscopic red- 
shifts, a number of ULIRGs are found in the ERCSC survey, but 
no HLIRGs, apart from the quasar 3C273. 

Fig. 13 shows the colour-ratio 5857/5545 versus redshift. 
Galaxies in the IRAS Bright Galaxy Sample are indicated as 
red dots. Sources for which T2 < 10 K in the 2-temperature 
fits discussed above are indicated as blue dots. Sources with 
logioiS ssj/S 54s) < 0.4 represent a new population of cooler 
submm sources. We have modelled all 17 galaxies with cool 857 
GHz/545 GHz colours (/ogio(5 857/5 545 < 0.4), and good opti- 
cal data (Fig. 14). Almost all require the very cold cirrus model 
with Tdust = 10-13 K,<p- 0.1 (green curves in Fig. 14). Table 1 
gives the properties of the galaxies whose SEDs we have mod- 
efled in detail. The columns are: ERCSC RA and Dec, IRAS ESC 
name, separation between ERCSC and IIFSCz position, spectro- 
scopic redshift, luminosity in cirrus (p - 5, (p - I and (p - O.l 
components, luminosity in M82 and Arp220 starburst compo- 
nents, optical luminosity, optical SED type, extinction, stellar 



9 



Planck Collaboration: The Planck View of Nearby Galaxies 



20 



o 



15 



10 








Fig. 10. Template fits for the four archetypal nearby galaxies, M51, MlOO, M82 and Arp 220. Black curves: fits with Efstathiou 
and Rowan-Robinson templates (black, separate components as dotted lines), blue curves: Silva et al. (1998) models. Planck 
ERCSC data shown as red filled hexagrams. /S'O-SWS mid-infrared spectroscopy data for M82 and Spitzer-IRS> data for Arp 220 
(Siebenmorgen & Kriigel 2007) are shown in magenta. 



mass and dust mass. Stellar mass and dust mass are calculated as 
in Rowan-Robinson et al. (2008), but there may be some model- 
dependencies in the dust masses since we have assumed that the 
cold dust in these galaxies is similar to Galactic dust. We can 
summarize the SED modelling shown in Figs 10, 11 and 14 as 
follows: (1) most nearby galaxies show evidence for dust at tem- 
peratures similar to that seen in the solar neighbourhood (0 = 1), 
as well as the warmer dust found in IRAS, ISO and Spitzer stud- 
ies; (2) there is a new population of cool submm galaxies with 
even cooler dust (</> = 0.1, Tf/,,,, = 10-13 K). This cooler dust is 
likely to have a more extended spatial distribution then generally 
assumed for the gas and dust in galaxies. 



6. Conclusions 

Our studies of nearby galaxies detected by Planck have shown 
evidence for colder dust than has previously been found in local 
galaxies. The temperature of this dust, however, is similar to that 



found in the solar neighbourhood. We also find that some local 
galaxies are both luminous and cool, with properties similar to 
those of the distant SMGs uncovered in deep submm surveys. 
This suggests that previous studies of dust in local galaxies have 
been biased away from such luminous cool objects. We also find 
that the dust SEDs in most galaxies are better described by para- 
metric models containing two dust components, one warm and 
one cold, with the cold component reaching temperatures as low 
as 10 K. Some objects have SEDs dominated by this cold mate- 
rial. These conclusions are based on both parametric fits and by 
detailed fitting of radiative transfer derived physical templates to 
the SEDs. However, other physical or parametric descriptions of 
dust, for example where /3 varies with wavelength, might lead 
to different results regarding this very cold component. This pa- 
per represents the first exploitation of Planck data for the study 
of a large sample of galaxies in the local Universe. As such it 
indicates both the benefits and hazards of the ERCSC for this 
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Fig. 11. Template fits for 8 nearby galaxies with detections in 16 bands. Blue curve is cirrus template used for solar neighbourhood 
(Tdust - 15-20 K, = 1), other components shown as dotted lines. 



work, but it also clearly demonstrates the tremendous potential 
of Planck data for the study of dust in galaxies. 
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Table 1. Parameters for Planck ERCSC-IIFSCz galaxies with SED fits 



RA 



Dec 



FSS dist z.pec, log(L) log(L) log(L) log(L) 

ciiT cirr crrr SB 

^ = 5 = 1 i/f = 0.1 



log(L) log(L) type 
Arp220 Opt 



Av log(M.) log(Mdi,s,) 



Template 



202.46930 47.19290 F13277+4727 0.1 0.001544 10.05 9.49 

M51 

185.72728 15.82777 F12203+1605 0.4 0.005240 10.53 10.19 

MlOO 

148.97540 69.68221 F09517+6954 0.2 0.000677 9.19 

M82 

233.74080 23.49795 F15327+2340 0.4 0.01813 

Arp220 



9.60 10.35 Sbc 0.0 10.75 7.54 

9.98 11.22 Sbc 0.0 11.62 8.15 

10.41 9.95 E 0.0 10.54 7.18 

12.15 10.58 E 0.0 11.17 8.15 



16-band, z <0.01 



175.24310 
143.03957 
213.89354 
181.50023 
206.59677 
184.95674 
218.18851 
178.20813 



11.47154 
21.50516 
36.22261 
47.47356 
43.86664 
29.61180 
49.45234 
44.12143 



Fl 1383+ 1144 

NGC3810 
F09293+2143 

NGC2903 
F14134+3627 

NGC5529 
Fl 2034+4745 

NGC4096 
F13443+4407 

NGC5297 
F12173+2953 

NGC4274 
F143 10+4940 

NGC5676 
Fl 1502+4423 

NGC3938 



0.1 



0.2 



0.3 



0.3 



0.5 



0.7 



0.4 



0.5 



0.003312 
0.001855 
0.009591 
0.001888 
0.008036 
0.003102 
0.007052 
0.002699 



9.70 
9.70 

8.90 
9.90 

10.60 



9.69 
9.49 
10.19 
9.09 
9.99 
9.64 

9.59 



9.24 
9.38 
9.60 
8.60 



9.10 



10.00 Sab 0.0 10.57 7.58 

10.41 Scd 0.2 10.79 7.43 

10.70 Sab 0.65 11.27 8.44 

9.65 Sab 0.6 10.22 6.95 

10.60 Scd 0.4 11.04 7.86 

10.40 Scd 0.0 10.94 7.45 

10.70 Sab 0.4 11.27 7.71 

9.90 Scd 0.5 10.34 7.40 



Cold log(S 857/5 545 )< 0.4 



159.78535 


41.69160 


F10361+4155 


2.0 


0.002228 






8.52 


8.10 


7.50 


sb 


0.2 


7.52 


7.33 


179.31494 


49.28748 


Fl 1547+4933 


0.4 


0.002592 




8.79 






8.70 


sb 


0.0 


8.72 


6.60 


218.19263 


9.88910 


F14302+1006 


0.6 


0.004574 


9.40 




8.92 




9.60 


Scd 


0.2 


10.04 


7.76 


208.22507 


-1.12087 


F13503-0052 


0.4 


0.004623 




9.64 






9.70 


Scd 


0.5 


10.14 


7.45 


148.79717 


9.27136 


F09521+0930 


0.0 


0.004854 






8.92 


9.40 


9.60 


Scd 


0.0 


10.04 


7.73 


219.79150 


5.25526 


F14366+0534 


0.4 


0.005020 




9.09 


8.62 


8.70 


10.35 


Scd 


0.2 


10.79 


7.54 


158.12950 


65.03790 


F 10290+65 17 


0.4 


0.005624 


9.40 




8.62 




9.80 


Scd 


0.2 


10.24 


7.48 


210.53130 


55.79348 


F14004+5603 


1.4 


0.006014 


9.05 




8.67 




10.25 


Scd 


0.7 


10.69 


7.50 


208.72266 


41.30995 


F13527+4133 


0.0 


0.007255 


9.51 


8.40 






9.96 


Scd 


0.2 


10.40 


6.76 


205.57553 


60.77630 


F13405+6101 


0.3 


0.007322 




9.20 


8.63 




9.61 


Sbc 


0.5 


10.01 


7.58 


228.37207 


58.49204 


F15122+5841 


1.0 


0.008474 


9.26 




8.95 




9.51 


Scd 


0.15 


9.95 


7.77 


231.67357 


40.55709 


F15248+4044 


0.4 


0.008743 




9.70 


9.43 


9.41 


10.16 


Scd 


0.7 


10.60 


8.31 


226.95305 


54.74415 


F15064+5456 


0.7 


0.01043 




10.00 






10.16 


Scd 


0.4 


10.60 


7.81 


126.57545 


22.88144 


F08233+2303 


1.3 


0.01794 




10.20 


10.13 




9.81 


Scd 


0.3 


10.25 


8.99 


231.43570 


52.44624 


F15243+5237 


0.3 


0.01948 


10.19 




9.81 




10.44 


Sab 


0.2 


11.01 


8.64 


237.68201 


55.60954 


F15495+5545 


0.2 


0.03974 


10.55 




10.52 


10.20 


10.70 


Sab 


0.0 


11.27 


9.34 


66.76954 


-49.12881 


F04257-4913 


0.8 


0.05828 






10.62 


11.55 


10.85 


Sab 


0.0 


11.42 


9.43 



Positions given are from the ERCSC. Where the sources have a well known name this is given beneath the IRAS naine. All luminosities are measured in solar 
luminosities. 
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Fig. 12. 545 GHz flux density versus redshift, showing loci for 
Arp220 template with Lir =10'^ and IO'-'Lq. Red dots are those 
galaxies in the IRAS Bright Galaxy Survey (BGS). We show 
here that the BGS does not sample as wide a range of galaxy 
properties as the Planck-IIFSCz sample discussed here. 
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Fig. 13. Sgs?/ S545 colour-ratio versus redshift. Points in red are 
IRAS BGS sources, points in blue are those sources identified as 
having T2 < 10 K on the basis of parametric fits. 
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= 0.1. 
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